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Abstract 

New results [] on the measurement of the proton structure function ^(x, Q 2 ) are reported 
for momentum transfers squared Q 2 > 1.5 GeV 2 and Bjorken x >3.5T0 -5 , using data collected 
by the HERA experiments HI and ZEUS in 1994. F-i increases significantly with decreasing x, 
even in the lowest reachable Q 2 region. The data are well described by a Next to Leading Order 
QCD fit, and support within the present precision that the rise at low x within this Q 2 range is 
generated via the DGLAP evolution equations. A comparison with models based on pomeron 
exchange is also presented. The gluon density is extracted and observed to rise at low x. 



1 Introduction 

The HERA ep collider has been designed to study Deep Inelastic Scattering (DIS) at very high Q 2 
where substructure of quarks could be observed. However in the first 3 years of data operation, 
which allowed a steady growth towards the design luminosity of the machine, most of the interest 
has focused on the study of low x, low Q 2 DIS, where new tests of perturbative QCD can be 
performed. The first observations on the 1992 data showed a rise of the proton structure function 
F2(x,Q 2 ) at low x < 1CU 2 with decreasing x || which was confirmed with the more precise 
data of 1993 |E|, ||]. Such a behaviour is qualitatively expected in the double leading log limit of 
Quantum Chromodynamics 0], It is, however, not clarified whether the linear QCD evolution 
equations, as the conventional DGLAP evolution || in InQ 2 and/or the BFKL evolution Q in 
ln(l/x), describe the rise of F2 or whether there is a significant effect due to nonlinear parton 



recombination [10J. At low Q 2 (< 5 GeV 2 ) the new results can be confronted to Regge inspired 
models, which expects a rather flat behaviour as a function of x, in order to study the transition 
between DIS and photoproduction. The 1994 data have allowed to reach Q 2 =1.5 GeV 2 and confirm 
the persistance of the rise at low x. The measurements have been achieved by using dedicated data 
samples (sect. 2) and are discussed and analyzed in terms of perturbative QCD in sect. 3. 



2 Structure Function Measurement 

In 1994 both experiments have reduced the minimum Q 2 at which they could measure F2 using 
several techniques: i) both experiments were able to diminish the region around the backward beam 
pipe in which the electron could not be measured reliably in 93, thus increasing the maximum polar 
angle of the scattered electron (measured w.r.t. to proton beam direction). This large statistic 

x Talk given at the 2 nd Rencontres du Vietnam, held in Ho Chi Min City in October 1995. The results presented 
here are extracted from the recent publications of the the HI and ZEUS collaborations 111 El. 



sample, taken with the nominal HERA conditions has an integrated luminosity of about 3 pb -1 . 
ii) An integrated luminosity of ~ 60 nb _1 of data was collected for which the interaction point was 
shifted by +62 cm, in the forward direction, resulting in an increase of the electron acceptance 
(so-called "shifted vertex" data sample ). hi) Both experiments used DIS events which underwent 
initial state photon radiation detected in an appropriate photon tagger to measure F 2 at lower Q 2 
(so called " radiative" sample) . The incident electron energy which participate in the hard scattering 
is thus reduced, and so is the Q 2 . The luminosity was determined from the measured cross section 
of the Bethe-Heitler reaction e~p — ► e~p r y, measuring the hard photon bremsstrahlung data only. 
The precision of the luminosity measurement is 1.5% (3.9% for the shifted vertex data). 

The kinematic variables of the inclusive scattering process ep — > eX can be reconstructed in 
different ways using measured quantities from the hadronic final state and from the scattered 
electron. The choice of the reconstruction method for Q 2 and y determines the size of systematic 
errors, acceptance and radiative corrections. The measurements presented here have been obtained 
with the electron (E) and with the E methods |0] for which the y, Q 2 and x formulae are 

E' e . 2 6 e E 2 E 2 sivL 2 9 e Qe,s 

Ve = 1 - ~et sin — y E = „ . F/n t-t Q e .s = -j ^e,s = 1) 

E e 2 E + E' e {l -cos0 e ) ■ l-y e ,E sy e ,s 

and E, p x ,Py,Pz are the four- momentum vector components of each particle, E e is the electron beam 
energy, s the squared center of mass energy of the collision, E = J2h ~ Pz,h and the summation 
is done over all hadronic final state particles neglecting their masses. The E method, which is 
independent of the hadronic final state, apart from the requirement that the interaction vertex is 
reconstructed using the final state hadrons, has at large y the best resolution in x and Q 2 but 
needs sizeable radiative corrections. At low y the E method is not applied due to the degradation 
of the y e resolution as 1/y. The E method, which has small radiative corrections, relies mostly on 
the hadronic measurement which has still an acceptable resolution at low y values and can be used 
from very low to large y values. HI measures F 2 with the E and the E method and after a complete 
consistency check, in particular at low x, uses the E method for y > 0.15 and the E method for 
y < 0.15. ZEUS measures F 2 at low Q 2 with the E method. 

The event selection is similar in the two experiments. Events are filtered on-line using calori- 
metric triggers which request an electromagnetic cluster of at least 5 GeV not vetoed by a trigger 
element signing a beam background event. Offline, further electron identification criteria are ap- 
plied (track-cluster link, shower shape and radius) and a minimum energy of 8(11) GeV is requested 
in ZEUS (HI). HI requests a reconstructed vertex within 3cr of the expected interaction position, 
while ZEUS requires that the quantity S = E + E' e (l — cos 9) satisfies 35 GeV < 5 < 65 GeV. If no 
particle escapes detection, 5 = 2E = 55 GeV, so the 5 cut reduces the photoproduction background 
and the size of the radiative corrections. The only significant background left after the selection 
comes from photoproduction in which a hadronic shower or a photon fakes an electron. In HI for 
instance, it has been estimated consistently both from the data and from Monte Carlo simulation 
and amounts to less than 3% except in a few bins where it can reach values up to 15%. It is 
subtracted statistically bin by bin and an error of 30% is assigned to it. 

The acceptance and the response of the detector has been studied and understood in great 
detail by the two experiments: more than two millions Monte Carlo DIS events were generated 
using DJANGO and different quark distribution parametrizations, corresponding to an inte- 
grated luminosity of approximately 20 pb _1 . The program is based on HERACLES for the 
electroweak interaction and on LEPTO [14| and ARIADNE [15] to simulate the hadronic final 
state. HERACLES includes first order radiative corrections, the simulation of real bremsstrahlung 
photons and the longitudinal structure function. For the parton densities, the GRV and the 
MRS parametrizations [|l7|] were used. The Monte Carlo events, after a detailed simulation based 
on the GEANT program, were subject to the same reconstruction and analysis chain as the real 
data. 
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Figure 1: Measurements of the proton structure function F2(x,Q 2 ) in the low Q 2 region by HI 
and ZEUS are shown together with the results from the E665 and NMC fixed target experiments. 
Different F2 parametrizations are confronted to the data: DOLA and CKMT (only in the upper 
row ofQ 2 bins); CTEQ3M, MRSG and MRSA' (lower row); GRV is shown for the full range. 

The structure function F2(x,Q 2 ) was derived after radiative corrections from the one-photon 
exchange cross section since effects due to Z boson exchange are smaller than 1% at low Q 2 . 

S " 2 ™\ { 2-2 v + T ^m X ,Q>) (2) 



dxdQ 2 Q A x v " 1 + R' 



The ratio R = F2/2xF\ — 1 was calculated using the QCD relation [18]. With the different 
data sets available detailed cross checks could be made in the kinematic regions of overlap. The 
results were found to be in very good agreement with each other for all kinematic reconstruction 
methods used, and the effect of systematic errors could be monitored: for the E method the main 
source of error are the energy calibration (known at the 1% level), the knowledge of the electron 
identification efficiency, the error on the polar angle of the scattered electron, (56 =lmrad) and 
the radiative corrections at low x. For the £ method, the knowledge of the absolute energy scale 
for the hadrons, the fraction of hadrons which stay undetected in particular at low x, due to 
calorimetric thresholds and to a lesser extent the electron energy calibration are the dominating 
factors. Further uncertainties common to all methods (selection, structure function dependence 
etc.) were also taken into account. The total F2 errors on the 1994 data ranges between 5 to 10% 
in the 10-100 GeV 2 range and between 10 to 20% below 10 GeV 2 . The final results on the 1994 
data of HI and ZEUS are shown in fig. [I] in the new kinematic domain reached using the radiative 
and the shifted vertex data. Compared to the 1993 data analyses the F2 measurement has been 
extended to lower x (from 1.8 • 10~ 4 to 3.5 • 10~ 5 ) and Q 2 (from 4.5 GeV 2 to 1.5 GeV 2 ). Both 
experiments are in good agreement and show that the F2 rise at low x persist, albeit less strongly, 
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down to the lowest measured Q 2 =1.5 GeV 2 . A smooth transition to the fixed target experiments 
E665 [19 1 and NMC [20| is observed with the low y results of HI, allowing to confront all these 
results to theoretical expectations. 



3 Low Q 2 and Perturbative QCD 



In fig. |] are also shown the extrapolations of the F2 parametrizations based on some theoretical 
model fitted to the previous data. They can be divided in two categories: one, motivated by Regge 
theory, assumes a pomeron exchange as a dynamical basis and successfully describes the behaviour 
of the total cross-sections of photoproduction and hadron-hadron collisions; the other is based on 
perturbative QCD and is known to well describe the DIS regime, but is expected to break down 
for a given x at some low Q 2 . The Regge models were expected to work at least at low Q 2 , but the 
DOLA parametrization which uses a "soft" pomeron (intercepts 1.08) [^] largely underestimate 
F2 at low x even at 1.5 GeV 2 . The CKMT model p3[ , which assumes that in the present Q 2 
range the "bare" pomeron becomes visible and has a higher trajectory intercept (~ 1-24), predicts 
a weaker rise at low x than observed, except maybe at 1.5 GeV 2 . These comparisons underline 
the difference between the behaviour of the total cross-section of real and virtual photons, since 



in the HERA kinematic domain and using the Hand 24] definition of the photon flux aj ^ can be 
expressed as a] a f(x,Q 2 ) ~ ^g^F 2 (x, Q 2 ). 

The parametrizations based on the DGLAP QCD evolution equations describe the data re- 
markably well, as expected above 5-10 GeV 2 , but surprizingly at values around 1 or 2 GeV 2 where 
non-perturbative effects were believed to distort the DGLAP picture. The MRSA' parametriza- 
tions of the parton densities are defined at Qq = 4 GeV 2 , then evolved in Q 2 and fitted to previous 
experimental data, including the 1993 HERA data. The agreement observed above 10 GeV 2 con- 
firms that the 1993 and 1994 HERA results are compatible. Between 1.5 and 10 GeV 2 the good 
description tells us that within the present precision perturbative QCD can be applied in this 
range. More striking is the confirmation of the pre-HERA prediction of the F2 rise at low x by the 




Figure 2: Variation of the exponent A from fits of the form F2 ~ x x at fixed Q 2 values. 
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GRV model [25] which conjectured at a very low energy scale (// 2 =0.34 GeV 2 ) that the proton is 
formed by valence-like partons as shown in fig. [| and that the DGLAP equations can be applied to 
generate "radiatively" the rise of the gluon and sea-quark density at low x, when evolving towards 
higher Q 2 . The HI and ZEUS results are very well described by the GRV model at low Q as can 
be seen in fig. || but also in the full HERA kinematic range, from 1.5 to 5000 GeV 2 p6[| . This 
success support the idea that the rise at low x is a direct consquence of the DGLAP equations, 
and that non-perturbative effects are relatively weak at low x and low Q 2 . The evolution with Q 2 
of the strength of the rise can be quantified by fitting an x~ x (or equivalently a W 2X , W being 
the invariant mass of the 7* — p system) function at fixed Q 2 to F2(x),x < 0.1. The values of A 
obtained by the fit in each Q 2 bin are displayed in fig ^ and clearly confirm the long time prediction 
made for asymptotic free field theories like QCD 0] of a rise of F2 at low x, and that the strength 
of this rise increases with Q 2 . With the present data, it is however not possible to know precisely 
this strength below 5 GeV 2 , thereby postponing a definite test of perturbative QCD in this region. 




1Q -4 W -2 W -4 1Q -2 1Q -4 1Q -2 W -4 jq-2 JQ-4 jq-2 



Figure 3: HI measurement (black circles) of the proton structure function F2(x,Q 2 ) as function 
of x in different bins of Q 2 . The inner error bar is the statistical error. The full error represents 
the statistical and systematic errors added in quadrature. The curve represent a NLO QCD fit to 
the HI, BCDMS (open squares) and NMC (open circles) data at Q 2 >5 GeV 2 . 



To make fully use of the new precision reached with the 1994 data, the HI collaboration has 
performed a Next-to-Leading Order (NLO) QCD fit on the HI, BCDMS and NMC data with 
the conditions Q 2 > 5 GeV 2 , and x < 0.5 if Q 2 < 15 GeV 2 to avoid higher-twists effects. The 
HI measurements which extend up to 5000 GeV 2 (they are shown up to 120 GeV 2 in fig. ||) 
were fitted successfully and allow to constrain the gluon density at low x. The parton densities 
were parametrized at Qq=5 GeV 2 , in particular the gluon was expressed with 3 parameters as 
xg(x) = A g x Bg {l — x) Cg . The quark and antiquark components of the sea were assumed to be 
equal, and u set equal to d. As determined in |^7|], the strange quark density was taken to be 
s = (u + d)/A. Further constraints were coming from the quark counting rules and the momentum 
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sum rules. For A the value of 263 MeV was taken |29|]. A detailed treatment of the F% errors 
propagation on the gluon density has been done, resulting in the error bands of fig |I]b which 
represent xg(x) at 5 and 20 GeV 2 . A variation of A by 65 MeV gives a change of 9% on the gluon 
density at 20 GeV 2 which has not been added to the error bands. The accuracy of this determination 
of xg is better by about a factor of two than the HI result based on the 1993 data [28]. 
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Figure 4: a) Parton densities (valence quarks (u v ,d v ), gluon (g) and sea quarks) of the GRV model 
at the initial energy scale fi 2 = 0.34 GeV 2 . b) Gluon density at 5 and 20 GeV 2 determined by a 
NLO fit to the H1,NMC and BCDMS data. The error bands represent the full error except for the 
uncertainty on A. 



A rise of the gluon density towards low x is observed which is related to the behaviour of 
F2 oc x~ x . Accordingly, the rise of xg towards low x increases with increasing Q 2 . Finally we 
can observe in fig. || that the data at Q 2 < 5 GeV 2 , which were excluded from the fit, are still 
well reproduced by the fit evolved backwards in Q 2 . More data at low x and Q 2 < 1 GeV 2 
are nevertheless needed to be able to test the hypothesis of a gluon density which would take 
the valence-like shape displayed in fig. ||a when Q 2 —>■ 0.3 GeV 2 , and more generally, to better 
understand the dynamics at low Q 2 and high parton densities. The HERA experiments, which 
have last year upgraded the capabilities of their backward detectors, will be able to reach these low 
Q 2 with the data taken in 1995 and 1996. 
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